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When a relativistically intense linearly polarized laser pulse is incident on an overdense plasma, a dense 
electron layer is formed on the plasma edge which relativistic motion results in high harmonic generation, 
ion acceleration and incoherent synchrotron emission of gamma-photons. Here we present a self-consistent 
analytical model that describes the edge motion and apply it to the problem of incoherent synchrotron 
emission by ultrarelativistic plasma electrons. The model takes into account both coherent radiation reaction 
from high harmonics and incoherent radiation reaction in the Landau-Lifshitz form. The analytical results 
are in agreement with 3D particle-in-cell simulations in a certain parameter region that corresponds to the 
relativistic electronic spring interaction regime. 


I. INTRODUCTION 

Gamma-rays and hard X-rays have become widely 
applied since their discovery. Most of present-day 
gamma-ray sources are based on the radioactive de¬ 
cay, bremsstrahlung, and backward Compton scattering; 
however, one needs to deal with radioactive materials or 
large-scale electron accelerators (linacs or synchrotrons) 
to use them and it limits their avaliability. The grow¬ 
ing demand for gamma-ray sources in numerous fields 
(medicine, radiography, nuclear physics) drives search for 
more convenient and accessible sources of hard X-rays 
and gamma-rays. In the recent decades an outstand¬ 
ing progress in laser technologies has been achieved, and 
now high-power (> 100 TW) short-pulse lasers have be¬ 
come commercially available. Due to ultrahigh electro¬ 
magnetic fields that they produce, lasers may be used 
to improve brightness and flux of gamma-ray sources. 
Namely, laser wake field acceleration (LWFA) is used to 
produce high-charge electron bunches for bremsstrahlung 
and Common source^^l® some other techniques are 
also proposecP^®^, one of them is incoherent synchrotron 
emission of plasmas lit by high-power laser pulses. 

When a relativistically intense laser pulse interacts 
with a target, the target electrons are expelled from 
atoms and accelerated to relativistic speeds, then the 
laser field forces them to emit photons due to Comp¬ 
ton scattering. If the laser field amplitude Eq is such 
that ao = eEQlmcuj ^ 1 (where c is the speed of 
light, uj is the laser field angular frequency, m and 
e > 0 are the electron mass and charge, respectively), 
the spectrum of emitted photons becomes synchrotron¬ 
like with a tail up to MeV or above. Earlier estima¬ 
tions and numerical simulations show that a large frac¬ 
tion of the laser pulse energy may be transformed into 
gamma-rays if laser intensity is high enough, e.g. per- 
cent^nHm 

for I > 10^^ W cm ^ and tens of percenfPHI^ 


Electronic mail: dmserebr@gmail.com 


for I > lO^^Wcm”^, so such gamma-ray sources look 
very promising (compared to linear Compton scattering- 
based sources with lower conversion efficiency, and to rel¬ 
atively complicated LWFA sources). But these phenom¬ 
ena are commonly studied by particle-in-cell codes, and 
currently a few analytical models are proposed. 

In this paper we present a self-consistent analytical 
model describing nonlinear electron motion in a thin 
dense layer arising at a laser-irradiated plasma edge in 
a specific parameter region. The corresponding laser- 
plasma interaction regime is referred as the relativistic 
electronic spring (RES), and is mostl y cons idered in the 
context of attosecond pulse generatiorP^^. In the RES 
regime a sizable portion of the laser energy is transferred 
into electron oscillationJI^ that seems preferable for effi¬ 
cient photon emission. The thickness of the electron layer 
in the RES regime may be as small as 1 nm, which makes 
coherent synchrotron emission feasible^^KUl xhe wave¬ 
length of 1 nm approximately corresponds to the pho¬ 
ton energy of 1 keV, therefore higher photon energies can 
be obtained only by the incoherent synchrotron emission 
(ISE). However, for laser intensities I > 10^^ W cm“^ 
average photon energies of single electron synchrotron 
spectrum are far beyond the limit of coherent emission; 
that’s why ISE can be very efficient. We compute radia¬ 
tion pattern and other properties of ISE in RES regime 
using the proposed analytical model for the plasma edge 
dynamics. 

Analytical results are verified by 3D particle-in-cell 
(PIC) simulations that take into account photon emission 
with radiation reaction and electron-positron pair pro¬ 
duction from hard photon^. Despite PIC simulations 
reveal more complicated electron dynamics than assumed 
in the model, the model describes well the plasma edge 
motion and the gamma-ray radiation pattern in a certain 
parameter region. Numerical simulations also allow us to 
distinguish various laser-plasma interaction regimes and 
compare gamma-ray generation efficiency between them. 
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II. ANALYTICAL MODEL FORMULATION 


To describe collective electron dynamics and hard pho¬ 
ton emission, let us introduce a self-consistent model of 
electron layer movement. Several major assumptions are 
pointed out. First, under the light pressure electrons of 
the irradiated plasma edge form a very thin layer (in 
comparison with the laser wavelength) which moves in 
such a way that laser pulse reflects from the target al¬ 
most completely, so we can neglect electron density per¬ 
turbations behind the layer. Second, we consider normal 
incidence of linearly polarized laser pulses and assume 
that the layer electrons move in the polarization plane. 
Third, only collective electron dynamics is considered (we 
suppose that the dispersion of individual electron char¬ 
acteristics inside the layer does not affect hard photon 
emission drastically and just results in a smoothing of 
the radiation pattern, photon spectrum, etc). The last 
assumption is that the ion motion is neglected: we re¬ 
strict ourselves to the case of few-cycle laser pulses so 
that the interaction time is short enough, or to the case 
of quite low laser intensity so that the ion displacement 
is negligible. 

The total force driving the layer electrons consists of 
the following parts: force caused by the incident laser 
field, force from self-generated electromagnetic fields (i.e. 
fields coherently emitted by the layer), force caused by 
electron-ion separation and radiation reaction force (i.e. 
force caused by ISE). Unlike the model of Ref. [T71 we 
do not adopt the requirement that the laser field behind 
the layer is completely compensated by the self-generated 
fields: this requirement make s it im possible to obtain 
Lorentz factor of the electronInstead, we start 
from the equations of motion that takes into account all 
the force parts mentioned above: 


dpx 

dt 

dpy 

dt 


— Ex "^yBz T Exx-) 

(1) 

— Ey -\- VxBz T Eyy^ 

( 2 ) 

E = Ez + Ep, 

( 3 ) 

B = B/ -1- Bp, 

( 4 ) 


where p is the electron momentum normalized to me, v 
is the electron velocity normalized to c, E and B are the 
electric and magnetic fields, respectively, normalized to 
mccj/e, indexes I and p denote laser and plasma fields 
(plasma fields are coherent emission and charge separa¬ 
tion fields), and is the radiation reaction force (caused 
by ISE). Here we assume that the laser pulse is incident 
along X axis and the polarization plane is the xy plane. 

The models similar to Eqs. 0-0 have been described 
in Refs. [20]and|2Tl They account only the radiation reac¬ 
tion caused by the fields coherently emitted by the layer, 
i.e., transverse components of Ep and Bp. On the con¬ 
trary, this model considers both the coherent part of ra¬ 
diation reaction (from coherent high harmonics emission) 
and the radiation reaction force caused by the incoher¬ 
ent emission of high-frequency photons. The wavelengths 



x/X 


FIG. 1. Schematic structure of the electron layer at the 
plasma edge and fields in RES regime. Curves: electron den¬ 
sity, gray; ion density, black dotted; Ex, green; sum of the 
laser field and the field emitted by the layer in -\-x direction 
{Ei^y Ey^-^), red; the wave emitted by the layer in —x direc¬ 
tion (Ey^-), blue. 


of these photons (emitted by single electrons) are much 
smaller than the layer thickness so they cannot be coher¬ 
ently summarized. 

As stated above, we consider the incidence of a laser 
pulse on a plasma half-space (unlike Ref. |2T| considering a 
foil much thinner than the laser wavelength). Therefore, 
the areal charge density of the layer is a function of the 
layer position, and according to RES assumptioiP^ is the 
following: 


e = noxi, 


( 5 ) 


where the charge density is normalized to crier '^o 
is the initial electron density normalized to the critical 
plasma density Ucr = mcj^/dTre^, xi is the layer displace¬ 
ment from the initial plasma edge position normalized 
to A/27r = c/cj. The fields coherently emitted by the 
layer can be then easily found from the Maxwell’s equa- 
tionP™ 


F - R - 
Ey-= 

’ 2{l + vC 

€X 

Ex = — \i < X < xi, 

X£ 


(6) 

( 7 ) 

(8) 


where indices + and — denote field values at ± 0 and 
also denote waves emitted by the layer in and —x 
directions. Ex can be found as a fiat capacitor field since 
the ions are immobile. An example of the fields generated 
by the electron layer oscillating according to Eqs. (10)- 
( p^ (see below) is shown in Eig. 

The radiation reaction force caused by ISE may be¬ 
come significant at field intensities of W/cm^ 
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so it shouldn’t be neglected for the sake of consis- 
tenc^l^^HSl main term in Landau-Lifshitz approx¬ 

imation of the radiation reaction force is the followin^^ 




47r Ve 


( 7 E + p X B)' - (p • E)' 


(9) 


where re = ^ jrxK? is the classical electron radius and 
7 = a/ 1 + is the electron Lorentz factor. The radia¬ 
tion reaction efficiently decelerates the emitting electrons 
and modifies their movement, however, we observe that 
the trajectory doesn’t change qualitatively. Anyway, the 
Lorentz factor and its distribution along the electron tra¬ 
jectory is affected by the radiation reaction that can yield 
modification of total power and radiation pattern of the 
emitted gamma-rays. 

It can be easily estimated from Eq. that for a simple 
circular trajectory in rotating electric field radiation re¬ 
action force is equal to the Lorentz force when the dimen¬ 
sionless field magnitude reaches athr = V^A/drrre ~ 400 
for optical wavelengths. For more complex trajectories 
radiation reaction force can be higher than for circular 
trajectory. Anyway, for more strict consideration ra¬ 
diation reaction force should be taken into account if 
ao > 100 - 400. 

Since a field affecting the layer can be found as a half¬ 
sum of fields on both layer sides, the system of equations 
that governs collective electron dynamics turns to be 


dpx 

dt 


noXj 
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dpy 

dt 


1 + 


V 
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— V. 

y^h 

+ Frx^ 


)^ly 

+ Fry^ 

dxi 


_ Px 

dt 
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_ Vy 

dt 

— '-'y — 
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( 10 ) 

( 11 ) 

( 12 ) 

(13) 


A. Some model properties 


In this subsection we neglect radiation reaction force 
F^ in order to demonstrate clearly a number of results of 
the model. First, the equation for Lorentz factor can be 
found from Eqs. (10)-(13) and is the following: 


^ = -{Ely + Ey,+ )Vy- 


noXjVx 

2(1 + 1 > 7 |) ■ 


(14) 


Thus, the electron Lorentz factor is determined by the 
field transmitted through the layer, Eiy + F^,+, and by 
the effective longitudinal field reduced by a factor of 1 + 
Vy^J. Lets assume that 7 ^ ^ 1 and \Eiy _^,+ | ^ uq, 
hence d^jdt dpxjdt. In this case Eqs. (10) and 


together with d{v‘^ + Vy)/dt ^ 0 yield: VyEiy + tiqXi 
Vx)l2 « 0. This equation is equivalent to the following: 

^ly E Ey^J^ ~ 0 . 
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FIG. 


2. (a) The layer position 

no = 400 (red curve) and 


given by the model for 
ao — riQ — 4UU (^red curvej and ao = no = 100 (blue 
curve). The dashed black curve corresponds to the layer mo¬ 
tion that yields a complete compensation of the incident field 
behind the layei^for no = ao- (b) The Lorentz factor of the 
above-mentioned trajectories. Thin grey curve corresponds to 
Eqs. (10)-(13) with Fr neglected, ao = no = 400. (c) Maxi¬ 


mal Lorentz factor gained by the layer electrons for ao = no at 
9 periods of the incident field vs the incident field magnitude 
according to the analytical model (blue curve with markers), 
and the same with neglected (thin grey curve), (d) The 
field reflected by the layer, ao = no = 400, at t = 9A/c. The 
gray area denotes a trajectory half period bounded by Vx ^ I 
points, vertical dark grey line denotes a time instant where 

'^X ~ 1 • 


Therefore, Eqs. (10)-(13) lead to the layer dynamics 
that almost completely compensates the incident field be¬ 
hind itself. This means that the trajectory obtained from 
the model is close to the layer trajectory obtained in the 
framework of Refs. im and El in which the layer trajec¬ 
tory is found from the strict requirement Eiy^Ey^j^ = 0 , 
that can be rewritten as followJ^^ 


dxi ~ 

dt ^^ly + 


However, the latter model assumes that Lorentz-factor is 
constant along the layer trajectory and doesn’t allow to 
find the radiation pattern since it highly depends on 7 / 
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FIG. 3. (a) The layer trajectory for ao = 240 and no = 320 obtained from Eqs. (10)-(13), color corresponds to 7 /. (b) 
Radiation pattern for the theory (solid) and PIC (dotted); ^ is the angle in the polarization plane, 'd = 0 corresponds to the 
direction of x axis, (c) Lorentz factor in the model 7 ^, and the average Lorentz factor ^pic of electrons with 7 > 0.05ao. 


distribution along the trajectory. 

Since the terms of Eq. 0 may be of the same order 
and the last term should not be neglected, Uq cannot be 
excluded from Eqs. ([Tq|)-(P^ by a normalization, and not 
the only parameter rio/ao governs the layer dynamics, 
but both ao and no- This means that the scaling law 
proposed in Refs. [T3l [TTl and \T7\ can be invalid for laser- 
plasma interactions in RES regime. 

The described model properties are illustrated in Fig.|^ 
where numerical solution of Eqs. (10)-([T^ is shown. 
First, it is seen that X£ given by the model of the present 
paper is close to the layer position given by models of 
Refs. [T 3 I and [TTl Second, due to the field compensation 7 ^ 
increases with time quite slowly and even for several laser 
periods 7 ^/ao remains < 1. Furthermore, 7 ^/ao drops 
with the instantaneous increase of ao and no- Third, ra¬ 
diation reaction doesn’t change the interaction picture 
considerably. 


B. ISE in the model 


Eqs. (10)-(13) don’t allow analytical solution and 


should be solved numerically. Solution depends on the 
parameters ao, no and laser pulse shape. 

The sample electron layer trajectory for typical mod¬ 
eling parameters is shown in Fig. |^a). The asymptotic 
trajectory resembles number 8 and is similar to a parti¬ 
cle trajectory in a linear-polarized EM wave (for a free 
particle it looks like number 8 only in a proper reference 
frame, though). For ao = 240 the maximum value of 
Lorentz factor is 120 and corresponds to the middle part 
of trajectory (where curve intersects itself); the minimum 
7 value is in the order of ten. The distribution of the layer 
energy over the trajectory yields the electron layer radi¬ 
ation pattern. Given that energy is ultrarelativistic over 
the whole period ( 7 ^ 1 ), radiation mechanism is almost 
purely synchrotron — practically all gamma quanta are 
emitted in the tangential direction. So we can use the fol¬ 
lowing formula for the emission power of a single electron 


at each instant of time: 

2e^c7^ 

where R is the curvature radius of the trajectory at the 
current point. For gamma quanta emitted by the whole 
layer we get 


Ilayer ^ 


(17) 


On-axis distribution of the hard photon density is shown 
in Fig. I^a), and the radiation pattern is shown in 
Fig. [|by On the curved part of the trajectory the 
gamma-rays are generated most efficiently due to small 
curvature radius, but the radiation direction changes 
rapidly at this point. It’s the central part of the trajec¬ 
tory that mostly contributes into the radiation pattern. 

The gamma ray radiation pattern shown in Fig |^b) 
corresponds to the trajectory from Fig. [^a). The dotted 
line is the pattern that was p roduced by PIC 3D numeri¬ 
cal simulations (see Sec. Ill for details). The pattern has 


two lobes in the plane of the laser E and k vectors. Lobes 
position depends on the laser pulse intensity; with greater 
intensity, the maximums are closer to the x-direction. In 
this model, nothing is emitted strictly forward because 
the layer electrons are never moving exactly in the -\-x 
direction. Qualitatively the model and numerical results 
are in agreement with each other, however simulations 
don’t give sharp peaks in the radiation pattern due to 
electron velocity spread which is neglected in model. 


III. RESULTS OF NUMERICAL SIMULATIONS 

In order to validate the model results, we have per¬ 
formed a series of numerical PIC simulations. The code 
used for is fully three-dimensional and can take into ac¬ 
count the quantum electrodynamical (QED) effect^. 
The code supports parallelization so it can be efficiently 
used on computational clusters. 
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FIG. 4. (a) On-axis electron and photon densities obtained 
in PIC for ao = 240 and no = 320. (b) For the same pa¬ 
rameters, the analytical trajectory (red) and xa according to 
Eq. |T^ for A = 0.01 A (lower boundary) and A = 0.11 A (up¬ 
per boundary), (c) Tracks of individual electrons obtained 
in PIC simulation. The color intensity corresponds to 
The average trajectory of the electrons with 7 > 0.05ao, xpic- 


the Xg distance from the pulse center. Therefore the con¬ 
figuration becomes much close to idealized modeling sit¬ 
uation (falling of a plain wave onto a flat target). 


A. Electron layer dynamics 


In a certain region of simulation parameters a layer 
with high electron density is formed in the front side of 
the target (see Sec. IV for details). Its thickness is in 
the order of 100 nm. Due to v x B force caused by lin¬ 
early polarized laser pulse the plasma edge is oscillating 
in the direction of pulse propagation (see Fig.|^a)). Elec¬ 
trostatic force caused by the ions compensates the light 
pressure on the average and holds the layer. 

From the numerical simulations we can acquire the 
layer trajectory and compare it to the analytical one. 
Current position xa of the layer in PIC is calculated us¬ 
ing the following formula 


A = 


f 


ne{x)lncr dx, 


(19) 


where A is some given in advance number of electrons. 
If the electrons form a thin layer and A is less than a 
number of electrons in the layer, xa is approximately 
equal to the layer position for any value of A. It is seen 
in Fig. I^b) that for 0.05 < A/A < 0.1 (shown by green), 
Xa obtained from Eq. ( [T^ is in good agreement with the 
analytical model. It is also seen that a small number of 
electrons can be pulled out from the plasma (see A < 
0.05A, shown in yellow, in Fig. |^b)). 

However, the layer dynamics in PIC simulations 
doesn’t fully obey theoretical assumptions. From 
Fig. S c) we can see that the layer is an average charac¬ 
teristic because single electrons constantly join and leave 
the layer. A typical electron lives in the layer one-two 
laser field periods and then enters the quasi-unperturbed 
plasma; it results in decrease of the layer energy com¬ 
pared to the theory (see jpic on Fig. [^c)) and modifi¬ 
cation of spatial distribution of energetic electrons (see 
xpic on Fig. Qc)). It alters the energetic properties 
of electrons and gamma-quanta in PIC compared to the 
theory. 


We have simulated the interaction of a falling 1 pm 
wavelength linear-polarized laser pulse with a fiat 2 pm 
wavelength foil. The laser pulse intensity and foil den¬ 
sity varied between different PIC runs. In this series, we 
considered the ions as fixed; conditions of applicability 
of this assumption are discussed in the section Laser 
pulse envelope has the following shape: 

= (18) 

It corresponds to a wavepacket which has almost constant 
amplitude in the central area and promptly decreases at 


B. Energetic properties and scaling laws 


One of well-known features of many laser-plasma prob¬ 
lems is that at relativistic laser intensities there is a sim¬ 
ilarity lavP^ the interaction properties depend on the 
dimensionless parameter 


S = 


Tlp_ 


CiQTlcr 


ao 


( 20 ) 


and not on no and uq separately. The current model 
doesn’t obey this scaling law strictly: equations ( TQp3 ) 
can’t be reduced to a form with one dimensionless pa¬ 
rameter even if radiation reaction is neglected (see also 


















6 



(^, rad 


that: first, electrons in PIC simulations are not monoen- 
ergetic and have a certain Lorentz-factor distribution, so 
there is always a number of electrons that are hotter than 
the layer electrons in the model. Due to the radiation law 
VP ^ 7 ^, their contribution into the total energy output 
can essentially increase the gamma-ray generation power. 
Second, the model estimation for the number of electrons 
in the layer Ue = noxi is obviously not true when layer 
disposition xi is about zero or even negative. In PIC 
simulations, the layer electrons sometimes had near-zero 
disposition, so model underestimates radiation of elec¬ 
trons in this part of trajectory. 


IV. CONDITIONS OF APPLICABILITY 


FIG. 5. Gamma-ray radiation pattern in PIG simulations at 
tto = n-o (aS = 1) varying from 40 to 480; lines show radiation 
pattern from the model for S' = 1, ao = 120 (solid line), 
ao = 240 (dashed). 


Sec. II A). However, the dependence on ao at S' = const is 
slight. We have compared radiation patterns from model 
and PIC simulations at different laser pulse intensities 
and foil densities but the same S parameters and the 
difference is negligible (see Fig. [^. 

Figure ^a) shows gamma-ray generation efficiency 
map in the no — ao axis which is acquired from 144 
PIC simulations; each point corresponds to the energy of 
gamma-rays generated during simulation timespan, nor¬ 
malized to the energy of the laser pulse. Different val¬ 
ues of similarity parameter S correspond to sloping lines 
starting at the coordinate origin, it can be clearly seen 
from the figure that S determines different regimes of 
gamma-ray generation. 

Higher ao values always result in increase of gamma- 
ray generation efficiency. The parameter region S < 1 
(which is denoted as ’T”) corresponds to the relativisti- 
cally self-induced transparency (RSIT) regime; the target 
becomes transparent because of the relativistic electron 
mass reduction and effective ao-times decrease of plasma 
frequency and increase of plasma density the electromag¬ 
netic waves can propagate througlP^. Penetrating laser 
pulse produces extreme heating of electrons and ions in 
plasma and causes emission of hard photons due to ISE. 
In this region gamma-rays are generated most efficiently; 
however, the model is not applicable here because of 
RSIT and absence of the single electron layer being con¬ 
sidered (we will discuss applicability conditions in more 
detail in the Sec. IV). 

For S' > 1, we can compare generation efficiency from 
PIC and model, the example of comparison under con¬ 
stant S = 1 and variable ao is shown in thej^b). From 
the model, the efficiency is about one order less that from 
PIC simulations. However, the scaling corresponds to a 
power law of the approximately same degree in both cases 
that turns to be an important result. Possible reasons for 


Based on 3D PIC simulations, we can distinguish a 
parameter region where the developed model can be ap¬ 
plied. The main assumption of the model is that the 
incident laser pulse pushes the plasma electrons which 
form a single thin layer that reflects the laser pulse well 
so electrons behind the layer are almost unperturbed. We 
perform the analysis of electron Lorentz-factor distribu¬ 
tion along longitudinal axis in different regimes (Fig. [^. 
It can be seen that the case of S' ^ 1 corresponds to the 
model best. For S < 1, RSIT leads to effective propa¬ 
gation of the laser pulse through the target, and several 
electronic structures with high values of Lorentz-factor 
are formed across the target. For S > 2, numerical sim¬ 
ulations show that electron dynamics significantly differs 
from the case of S ^ 1 (where individual electrons stay 
in the layer up to several laser periods): the electrons 
generally escape the layer on each of laser half-periods, 
and new (‘cold’) electrons from the plasma replace them. 
Therefore electrons are grouped in thin bunches which 
propagate in x-direction (see Fig.[^ (c)). This effect leads 
to effective increase of overall electrons energy and de¬ 
crease of gamma-rays generation efficiency because elec¬ 
trons in the layer do not reach high energies during laser 
half-period. 

Therefore we define the region where the model can be 
applied as 1 < S' < 2. Limitation on ao parameter can 
be found from the requirement that we neglect quantum 
electrodynamical effects in this model so ao shouldn’t be 
greater than ^ 500. Obviously the laser pulse should be 
also strong enough so that the electron layer with high 
Lorentz-factor can be formed at all (thus ao ^ 1). 

One should consider attainable plasma densities as 
well. For solid targets a typical electron density is in 
the order of 10^^ — lO^^crn”^ or higher so no is in the 
order of hundreds (if A 1pm). In order to achieve lower 
densities, some exotic targets are required. 

The model also has the requirement of immobile ions. 
The numerical simulations that have been presented 
above are carried with immobile ions too, but sev¬ 
eral numerical experiments with moving ions were done. 
Charge-to-mass ratio of the ions was 0.25 of that for hy¬ 
drogen atoms. The results (gamma-ray radiation pat- 
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FIG. 6. (a) Gamma-ray generation efficiency map obtained 

from 144 PIG simulations. Foil thickness is 2 pm, foil density 
and laser pulse intensity are variable parameters. II is the 
region where the model is applicable; III is the region of rela- 
tivistically self-induced transparency (RSIT), in region I the 
stable electron layer is not formed in PIG. (b) Gamma-ray 
generation efficiency from PIG and analytical model at S' = 1 
(ao = no). 



XXX 


FIG. 7. Electron Lorentz-factor distribution along x-axis at 
t = 5.8A/c (laser pulse is incident from the left) for different 
values of S. Only electrons from the region near the x axis, 
where laser pulse front is plane, are shown. 


tern, electron layer dynamics) didn’t change drastically 
in case of moving ions, although gamma-ray generation 
efficiency became lower about 2 times. To get mostly 
aligned with the model and improve gamma-ray genera¬ 
tion efficiency, one should take as smallest ion charge-to- 
mass ratio as possible. 


V. CONCLUSION 

The self-consistent model that describes electron dy¬ 
namics in the interaction of a relativistically strong 
laser field with over dense plasma halfspace is presented. 
Plasma electrons are pushed by v x B force and form 
a thin layer on the plasma edge. This layer coherently 
emit electromagnetic spikes as the laser pulse reflection 
and considerably suppress the fields transmitted into the 
plasma. The back-reaction of the coherently emitted 
fields is taken into account, as well as radiation reaction 
caused by the incoherent emission of the layer. 

The back-reflected field presents a train of electromag¬ 
netic spikes emitted when the electron layer has maximal 
velocity towards the laser pulse. These time instances 
correspond to the minimal value of the electron Lorentz 
factor that means not optimal conditions of short-spike 
generation and low value of the cut-off in the spectrum 
of the reflected fielcP^^^. 

When electrons are ultrarelativistic, the supposition 
that the plasma dynamics is governed by the only pa- 
ram eter S = no /ao (not by no and ao separately) is often 
^sejnEEll. It yields the electron Lorentz factor 7 oc ao- 
But the presented theory demonstrates that the process 
of electron acceleration can be a fine effect. Namely, it 
occurs that dy/dt = 0 if we neglect the electron rest mass 
and fields behind the layer in the framework of the pre¬ 
sented theor}^^. The accurate analysis shows that the 
equation for the electron energy depends on two param¬ 
eters no and ao separately, and yields that the ratio 7 /ao 
drops with the increase of ao- The losses caused by the 
incoherent photon emission strengthen this drop. 

The presented model allows one to find electron trajec¬ 
tories (including Lorentz factors) and, hence, incoherent 
synchrotron emission from the plasma edge. Plasma edge 
dynamics and gamma-ray radiation pattern are in a fairly 
good agreement with the results of 3D PIC simulations. 
Namely, the model predicts a two-lobe radiation pattern 
with lobes approximately parallel to the electric field di¬ 
rection of the incident laser pulse. This radiation pat¬ 
tern is caused by the 8-like electron trajectories extended 
along this direction. It is also predicted that for S' « 1, 
in contrast to the electron reinjection regim^I^, emission 
towards the incident laser pulse is much weaker than in 
perpendicular direction, that is also in good agreement 
with 3D PIC simulations. Therefore, angular distribu¬ 
tion of gamma-rays can specify the laser-plasma interac¬ 
tion regime in experiments at ultrahigh intensity. 

PIC simulations reveal that a thin electron layer is 
formed on the plasma edge if S = no/ao ~ 1. In the case 

























of S' < 1/2 relativistically self-induced transparency oc¬ 
curs, and the laser pulse propagates through the plasma. 
If S > 2, the skin depth (the thickness of the current 
distribution in the plasma) becomes larger than electron 
displacement perpendicular to the plasma surface, and a 
lot of electrons escapes to the plasma bulk. If the laser 
pulse duration or the laser pulse magnitude uq is high, ion 
motion becomes considerable. This case is also beyond 
the region where the model is applicable. 

From numerical simulations we also see that individ¬ 
ual electrons stay near the plasma edge generally a half 
of the laser period, thus, electrons permanently leave and 
join the layer. Electrons that leave the layer transfer the 
laser energy into the plasma bulk and heat the target. 
Electrons that join the layer have energy which is sig¬ 
nificantly different from the average electron energy in 
the layer that may cause a sizable discrepancy between 
model predictions and PIC results. Namely, PIC and the 
model gives different values for overall electron energy 
and generation efficiency of gamma rays. Nevertheless, 
the analytical model yields the same scaling of the gen¬ 
eration efficiency as PIC if S' ~ I. 

Therefore, the presented theoretical model can be used 
for the analysis of coherent and incoherent photon emis¬ 
sion for no ^ no- The model can be used for optimization 
of high harmonic and gamma-ray generation in the in¬ 
teraction of ultraintense laser pulses with plasmas. Some 
results of this work would be extended to the case of 
obliquely incident laser pulses. 
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